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ABSTRACT 



Context. Clumpy disk galaxies in the distant universe, at redshift of z > 1, have been observed to host several giant clumps in their 
disks. They are thought to correspond to early formative stages of disk galaxies. On the other hand, halo objects, such as old globular 
clusters and halo stars, are likely to consist of the oldest stars in a galaxy (age > 10 Gyr), clumpy disk formation can thus be presumed 
to take place in a pre-existing halo system. 

Aims. Giant clumps orbit in the same direction in a premature disk and are so massive that they may be expected to interact gravita- 
tionally with halo objects and exercise influence on the kinematic state of the halo. Accordingly, I scrutinize the possibility that the 
clumps leave a kinematic imprint of the clumpy disk formation on a halo system. 

Methods. I perform a restricted A'-body calculation with a toy model to study the kinematic influence on a halo by orbital motions of 
clumps and the dependence of the results on masses (mass loss), number, and orbital radii of the clumps. 

Results. I show that halo objects can catch clump motions and acquire disky rotation in a dynamical friction time scale of the clumps, 
~ 0.5 Gyr. The influence of clumps is limited within a region around the disk, while the halo system shows vertical gradients of 
net rotation velocity and orbital eccentricity. The significance of the kinematic influence strongly depends on the clump masses; the 
lower limit of postulated clump mass would be ~ 5 x 10 s M Q . The result also depends on whether the clumps are subjected to rapid 
mass loss or not, which is an open question under debate in recent studies. The existence of such massive clumps is not unrealistic. I 
therefore suggest that the imprints of past clumpy disk formation could remain in current galactic halos. 

Key words, methods: numerical - galaxies: halos - galaxies: spiral - galaxies: evolution 
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Although disk galaxies in the local universe generally 
have relatively smooth stellar distributions in which small 
star clusters of M c i < 10 3 M Q reside (e.g., lLada & Ladal 
I2003I) . some galaxies in the high redshift universe (red- 
shift of z ~ 1 - 5) have been observed to have 
clum py structures in which massive clumps are f ormin 



(Nog uchl 119981 Il999fc llmmeli et al.l I2004allb1: iBournaud et all 
2007; Elmegreen et al. 2008a b; Bournaud & Elmegreen 2009; 



Agertz et all 120091: iDeke 



Cever ino et alJl2010L I20T2T 



et alJ l2009at jAumer et al 
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(e.g., van den Bergh et al. 1996; Elmegreen et al. 2004, 200 
2009a b; Elmegreen & Elmegreen 2005; Forster Schreiber et al 
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2006, 2009, 2011; Genz el et alj|20 06. 2008, 2011; Tacconi etail 
| 20ialPtiechll2010l:lGuo et al]l2012HSwinbank et alJl201 lLl2012t 
iNewman et alJl2012i: IWuvts et al.ll2012t iLivermore et aljbom 
Regardless of their irregular morphologies, the galaxies in- 
dicate clear rotation in spectroscopic observations as signa- 
tures of disk structures, although some c onsiderable fraction 
of them may be ongoing mergers (e.g., IWeiner et al.l [2006; 
iForster Schreiber et alJ 120091 : lPuechll2bl0h . These are referred 
to as clump clusters and chain galaxies, depending on viewing 
angles seen face-on or edge-on, respectively. In such galaxies, 
a single clump was observed to hav e, at the largest, a mass 
of M c i ~ 10 9 - 10 10 M . Recently, Ho dge et alJ (l2012l) also 
found a rotating clumpy gas disk in a sub-millimeter galaxy 
at a redshift z - 4. The clumpy galaxies are thought to be 
disk galaxies in early stages of their formation. Numerical sim- 
ulations showed that dynamical instability due to the gas-rich 
nature of a premature disk can lead to clump formation: the 
clumps fall into the galactic center by dynamical friction and 
may form a central bulge. Such numerical studies demonstrated 
that a clumpy galaxy can finally evolve into a disk galaxy 



However, these results do not necessarily mean that all 
of the current disk galaxies were once clumpy galaxies. 
Elmegreen & Elmegreenl (120051) have observationally discussed 
that clumpy galaxies have too high a surface density to be current 
disk galaxies such as the Milky Way and that they would evolve 
into other types of galaxy (e.g., ell iptical galaxies by merg- 
ers). Furthermore. iMartig et al.l d2012l) performed simulations of 
modeled galaxies joined to a cosmological run and demonstrated 
that clumpy disk formation is preceded by intense gas accretion 
at high redshift and that clumpy disk formation would not be 
common in isolated environments. 

How common was clumpy disk formation in the universe? 
Can we know which disk galaxy was once a clumpy galaxy? To 
answer these questions, it is important to seek a clue concern- 
ing the past clumpy disk formation in current galactic structures. 
In this paper, I focus on halo objects such as old globular clus- 
ters (GCs) and halo stars (HSs), which are likely to consist of 
the oldest stars in a galaxy. Generally, GCs are classified into 
metal-poor and -rich ones. Metal-poor GCs are usually old (age 
> 10 Gyr), although there are some exceptions and sub-groups 
in their classification (e g.. [Harris & Racindl 19791: iHarrii ll 991. 
l2001l:lBrodie&Stra der 2006, and references therein). 

Clumpy disk galaxies, on the other hand, can be observed 
even at a redshift of z ~ 1 (~ 8 Gyr in look-back time). Hence, 
clumpy disk formation can be thought to postdate the forma- 
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tion epoch of halo objects and expected to take place in a pre- 
existing halo if the halo objects have an internal origin, such as 
the collapse of pro to-galactic gafl which was the concept of 
Egge n et all {T962). If clumps are many and massive enough 
to interact gravitationally with halo objects, orbital motions of 
the clumps may exercise influence on th e ki nematic state of 
the ha lo objects. Elmeg reen et alJ (l2008ah and llnoue & S aitoh 
d201lh have conducted numerical simulations that showed that a 
density profile of a dark matter (DM) halo can be changed due to 
kinematic heating by giant clumps. Their results also imply that 
kinematics of halo objects may be affected by the clumps. 

In this paper, I perform restricted A^-body calculations us- 
ing simple toy models, demonstrate how giant clumps exercise 
influence on halo objects, and examine dependence on masses, 
number, and orbital radii of the clumps. Additionally, I discuss 
the impact of mass loss of the clumps on my result, which have 
been recently suggested to be caused by a strong gas outflow. 
Although, in my study I discuss the kinematics of GC systems, 
this can also be applied to HS. I explain my model settings and 
method in ^2]and show the results and discuss the kinematic im- 
print of clumps in |3] I present the discussion and summary of 
this paper in $4] and $5] respectively. 



2. Restricted A/-body calculation 

As mentioned above, I perform restricted A^-body calculations 
with a simple model. I assume a DM halo r epresented by a rigid 
potential of a Navarro-Frenk- White model dNavarro et alJ l 1997. 
hereafter NFW) with a virial mass of 1.0 x 10 12 M Q , a virial 
radius of 200 kpc, and a concentration parameter of 15. Indeed, 
clumpy galaxies at a redshift of z ~ 2 are thought to live in DM 
halos of ~ 10 12 Mq in which cold stream is expected to take 
place (e.g. JDekel etai1l2009bh . 

Initially GCs (and HSs) are assumed to have a distribution 
with a spherical number density profile of 



«Gc(r) = no 



ro 



3.5 



r 2 (r + r) 



1.5 ' 



(1) 



accor ding to observations of GCs in the Milky Way dHarrisI 
120011) . where I set ro = 4 kpc and no is a constant normalized 
by the total number of GCs (see below). I assume a truncation 
on the distribution outside 100 kpc and a constant density inside 
100 pc. The GC system has an equilibrial and isotropic kine- 
matic state determined by a solution of Jeans equations. A single 
GC is represented by a single mass-less particle. 

As a fiducial model, the structure of a clump is represented 
by a Plummer model with a m ass of M c \ = 5 x 10 8 Mr and a 
Plummer radius of e = 200 pc. Elmegreen & Elmegreen (2005) 
have observed that the numbers of clumps are 5 - 14 and the 
averaged clump masses range from 0.21 - 1.36 xlO 9 M in 
their sample galaxies. In accordance with their observation, my 
initial condition sets ten clumps in circular orbits at radii of 
1, 2, 3, • • ■ , 10 kpc from the halo center. All clumps orbit in the 
same direction in the same orbital plane. The initial orbital phase 
is randomly determined on each clump. Additionally, I perform 
other runs with different settings of the clump mass, number, and 
the initial radii (see 



1 As another scenario, halo objects may be brou ght in by accretio n 
events of dwarf galaxies after the clum py phase (e .g., Abadi et al. 2006). 
Composite models are also possible. Mac kev & Gilmora 12004) have 
estimated that nearly half of the mass of the Galactic stellar halo is 
from external accretion and the other half is from an internal origin. 




Fig. 1. Time evolution of orbital radii of the clumps in the fidu- 
cial model. 



Gravity calculation takes into account gravitational forces 
from the clumps to the GCs and mutual interactions between 
the clumps. I set the total number of GCs to 100,000. Of course, 
there is not such a large number of GCs in any given galaxy. 
However, results of this calculation are independent of the num- 
ber of GCs since they are mass-less particles. Hence, I employ as 
large a number of GCs as possible for the sake of good statistics. 
Dynamical friction is ignored on the GCs, but embodied with 
the Chandrasekhar formula on the clumps, adop ting In A = 8.4 
(IChan drasekhar 1943l: lBinnv & Tr emaine 2008). Sin ce a proper 
value of In A is unclear (e.g., lHashimoto et al.ll2.003l: [Fuiii et all 
2006: llnouel2009ll201 ll) . I set the value to match infall time scale 
of the clumps w ith the results of other numerical simulations , 
~ 0.5 Gvr (e.g.. | Bournaud et al.ll2007l: lElm egreen et al. 2008a; 
ICeverino et al.lf2012h . When a distance between two clumps be- 
comes shorter than e, they merge into one, the mass of which 
is the sum of the two. I conduct the orbital calculation with the 
second-order leap-frog time integrator and a shared time step of 
At = 0.03 x Ve 3 /(GM cl ) = 0.057 Myr until t = 5 Gyr. 



3. Results 

3. 1 . The fiducial model 

In the fiducial model, I set ten clumps in the DM halo; the masses 
of the clumps are M c \ — 5 x 10 8 M H , following the observation 
of Elmegreen & Elmegreen (2005). Since their galaxy samples 
have disk diameters of 13 - 29 kpc, I set the initial clump orbital 
radii to 1 - 10 kpc with constant separations. I assume that the 
clump masses are constant unless they merge. 

The clumps fall into the galactic center by dynamical friction 
and finally form a clump-origin bulge. FigureQ]shows time evo- 
lution of orbital radii of the clumps in the fiducial case. Nine 
of the ten clumps reach the galactic center by t 0.4 Gyr 
while merging with one another. Although the outermost clump 
does not experience any merger, it spirals into the center at 
t^ 1.2 Gyr. 

Figure [2] displays a GC distribution seen edge-on in the end 
state of the calculation. The figure shows a slightly aspheri- 
cal distribution, though it was initially assumed to be spherical. 
Figure[3]shows a line-of-sight (LoS) velocity map. The halo sys- 
tem indicates a clear rotation around the disk plane in the end 
state, implying that the GCs catch the clump orbital motions 
through gravitational interaction. 
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Fig. 2. Probability density distribution of the GCs from the edge- 
on view in the unit of kpc 2 in the end state of the fiducial cal- 
culation. The initial distribution was set to the spherical one, ac- 
cording to Eq. Q] The probability density is defined as the num- 
ber of GCs in each bin divided by the total number. 
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Fig. 3. Map of mean LoS velocity from the edge-on view in the 
fiducial run. The values are normalized by the LoS velocity dis- 
persion inside ro = 4 kpc, which is 106 km s _1 . 



Vertical profiles of mean azimuthal velocities and orbital ec- 
centricities in different radial ranges are shown in Figure [4] The 
kinematic influence by the clumps can be seen up to R < 10 kpc 
and z < 3 kpc, showing vertical gradients of the rotation veloc- 
ities and the orbital eccentricities. In this region, the halo indi- 
cates slow but clear net rotation and GCs rotating prograde are 
circularized. In the innermost region, R = 1 - 4 kpc, orbits of ret- 
rograde GCs are radialized rather than circularized, and prograde 
GCs are not strongly circularized. This is probably because of a 
potential of the massive bulge built up from the clumps. 

In the fiducial run, I see that clump motions can rotate the 
halo objects if they are many and massive enough. However, 
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Fig. 4. Mean azimuthal velocity (top) and orbital eccentricity, 
e = (r apo - r peri )/(r apo + r peri ), (bottom) of the GCs as functions 
of distance from the clump orbital plane in the end state of the 
fiducial run, where r apo and r pe ,„ are apo- and peri-center dis- 
tances, respectively. Each line indicates a radial range in cylin- 
drical coordinate, R. In the bottom panel, GCs rotating prograde 
and retrograde are separately plotted; the thick and thin lines 
correspond to prograde and retrograde ones, respectively. 



some parameters assumed in the calculation are unclear in cur- 
rent observations and clumpy galaxies have been observed to 
have a variety of their properties. In the following subsections, I 
perform further calculations with different settings. 

3.2. Fewer clumps 

Although Elmegreen & Elmegree^ (I2005I) have observed that 
the average number of clumps is nearly ten in a clumpy 
galaxy, this may not be a typi cal number. For example, 
iForster Schreiber et alj (|2011|) and iGuo et alj d2012l) observed 
that the median numbers of clumps are five and four in their sam- 
ples, respectively. Accordingly, I run similar calculations with 
fewer clumps. 

I set five clumps in the halo model, of which initial orbital 
radii are 2, 4, • ■ • , 10 kpc. The other parameters are the same as 
the fiducial run. Figure [5] displays a LoS velocity map and ver- 
tical profiles of mean azimuthal velocities and orbital eccentric- 
ities. Although the rotating signature is a little weaker than the 
fiducial run, significant differences cannot be seen. Hence, the 
number of clumps does not seem to change the results in the 
range of 5 - 10 clumps. 

The orbital radii of clumps in the fiducial run may also not 
be typical. Although giant clumps seem to fall t oward the cen- 
ter from > 10 kpc in co s mol ogical simulations dCeverino et al.l 
1201 2l) . lElmegreen et al.l (120071) had observed that root mean 
squares of clump positions in galaxies are 1 - 8 kpc. Thus, some 
clumpy galaxies may be more compact than the fiducial model. 
Accordingly, I conduct a run setting five clumps at 1 , 2, • • • , 5 kpc 
from the galactic center. 

The results are shown in Figure [6] The figure indicates sig- 
nificant rotation and orbital circularization in R = 1 - 7 kpc as a 
kinematic imprint of the clump motions. The region ofR>l kpc 
remains intact because of the absence of outer clumps. From 
these results, I suggest that inner regions of a disk galaxy may 
retain the disky rotation caused by clump motions since clumps 
have to pass through the inner regions of galaxies to reach the 
center, the inner regions therefore cannot avoid the kinematic in- 
fluence by the clumps, in cases where the clumps are long lived 
and can reach the galactic center by dynamical friction. 
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Fig. 5. Same as Figure [3] and Figure [4] but setting five clumps Fig. 6. Same as Figure [3] and Figure [4] but setting five clumps 
with constant separations of 2 kpc in the initial state. The LoS with constant separations of 1 kpc between R — 1 - 5 kpc in the 



velocity dispersion inside ro is 92.2 km s in the top panel. 



initial state. The LoS velocity dispersion inside ro is 92.5 km s 1 
in the top panel. 



3.3. Lower mass clumps 

Masses of giant clu mps have been observed to range widely in 
10 5 - 10 10 Mm (e.g..lElmegreen et al]|2009albt lGuQ eTaIll20Tl 
iLivermore et al.l2012l) . ISwinbank et all (120121) have recently dis- 
cussed that the most massive clump mass in a clumpy galaxy 
would be 10 s - 10 9 M and present observations have shown 
that massive clumps reaching > 10 9 M do exist. However, 
such massive clumps may be rare. Accordingly, I run another 
simulation setting ten clumps, the masses of which are M c \ = 
1 x 10 8 M . The other parameters are the same as the fiducial 
run. 

The result is shown in Figure [7] In comparison to the previ- 
ous runs, the halo rotation due to the clumps is much weaker in 
this case. In the bottom panel, although a rotating signature can 
be seen weakly in z < 1 kpc and R = 1 - 7 kpc, it is less sig- 
nificant than the previous cases and the orbital eccentricities are 
unchanged. Hence, clumps smaller than 1 x 10 8 M could not be 
expected to leave the kinematic imprint on halo objects. 

It is important to remark that dynamical friction depends 
on clump mass. In the Chandrasekhar formula, the decelera- 
tion is proportional to the clump mass an d the low-mass clumps 
slowly migrate into the galactic center JChandrase kharl [T943t 
Binnv & Tremaine 2008). Since the clump masses are five times 
smaller, the time scale of orbital shrinkage of the clumps is 
approximately five times longer than that in the fiducial case. 



Although such lower mass clumps have a longer time to interact 
with the halo objects, the imprinted rotation on the halo is much 
weaker. Accordingly, I see that the clump interaction with halo 
objects seems to strongly depend on the clump masses. 

3.4. The case with mass loss 

lNoguchil (il996) has pointed out the possibility that giant clumps 
may be fragile objects if they suffer mas s loss from outflow 
stron g enough to b l ow ga s out. Recently, iGenzel et al.l (1201 lh 
and iNewman et al.1 (120121) have actually observed such out- 
flow from clumps in some galaxies, which exceeds escap e ve- 
locities of the clum ps. Furthermore, Hopk ins et al.l (1201 ll) and 
Genel et 

ID 

d2012l) performed numerical simulations that sug- 
gest that such superwinds can blow gas away from clumps 
and prevent the clumps from growing massive. As a result, 
the clumps become vulnerable to tidal disruption. These stud- 
ies discussed that clumps would be short lived and disrupted 
in a f ew dynamical times of the c lump ro tation. On the other 
hand, iForster Schreiber et all d201ll) and iGuo et alJ d2012h ob- 
served radial-age gradients of clumps in some galaxies. They 
argued that the age gradients could be indicativ e of the migra- 
tion an d longevity of the clumps. Moreover, Elm egreen et al.l 
(2009a) discussed that clump ages are intrinsically closer to ages 
of bulge-like objects in galaxies, in cases where there are bulge- 
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Fig. 7. Same as Figure [3] and Figure |4] but ten clumps of 1 x Fig. 8. Same as Figure [3] and Figure [4] but the mass loss from 



10 M . The LoS velocity dispersion inside ro is 81.2 km s 



the clumps is implemented. The values are normalized by the 
velocity dispersion inside ro = 4 kpc, which is 75.0 km s _1 . 



like objects. They discussed that giant clumps are long lived and 
can reach the galactic center by dynamical friction. 

In cases where the superwinds are real and clumps are 
short lived, clumps may not be able to interact with halo ob- 
jects. Accordingly, I perform a run assuming mass loss from the 
clumps in the calculation to study how the mass loss caused by 
the superwinds and/or tidal disruption affects the kinematic im- 
print on a halo. In this run, the masses of the clumps are initially 
M c i = 5 x 10 8 M . The clumps decrease their masses at a rate of 
0.743 % per time step, while their orbits are time integrated. At 
this mass loss rate, the clumps lose 99 % of their masses in the 
first t = 0.35 Gyr (one orbital period at 10 kpc from the center)Q 
In this case, the clumps do not build up a bulge because of the 
mass loss. The initial condition have ten clumps and the other 
settings are the same as the fiducial model. 

Figure [8] displays a LoS velocity map, mean azimuthal ve- 
locities, and eccentricities in this case. Although halo rotation 
appears only weakly inside R < 4 kpc and z < 2 kpc, the weak 
influence in the innermost region may be due to the initial set- 
up arbitrarily putting some clumps on the inner radii. Here, the 
mass loss time scale becomes relatively long compared to orbital 
time scale. As noted in 33.31 since the deceleration by dynamical 
friction is proportional to the clump mass in the Chandrasekhar 
formula, the clumps losing their masses can keep their orbits for 
a long time. However, as shown here, the clumps subjected to the 

2 This rate corresponds to -6.56 M Q yr~' for M A = 5 x 10 s M Q . 



rapid mass loss do not seem to exercise gravity on the halo ob- 
jects significantly. Thus, I suggest that the existence of the rapid 
mass loss from clumps could be a crucial factor of the kinematic 
imprint of clumpy disk formation, which is still under debate. 

4. Discussion 

From the results above, I suggest that if giant clumps in a clumpy 
disk galaxy can retain their masses of > 5 x 10 8 M until reach- 
ing the galactic center, they could gravitationally interact with 
halo objects, such as old GCs and HSs, even in the time scale of 
dynamical friction, t 0.5 Gyr. Hence, clumpy disk formation 
can be expected to leave a kinematic imprint on the halo objects, 
which may still remain in a current galactic halo if the galaxy 
has not experienced violent disturbances, such as major merg- 
ers after the clumpy phase. In such a case, the influenced halos 
would show net rotation with vertical gradients of mean velocity 
and orbital eccentricity around a disk plane. 

The calculations, however, seem to show that the signifi- 
cance of the kinematic influence strongly depends on the clump 
mass. If the clumps are less massive than 1 x 10 8 M or if they 
suffer violent mass loss, the kinematic state of the halo barely 
changes. This means that clamp mass (mass loss rate) is a key 
to whether or not we can deduce remnants of clumpy galaxies 
from halo kinematics in current disk galaxies. From the result 
of 33.31 1 estimate that the critical mass to exercise influence on 
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a halo would be M c ; ~ 5 x 10 M Q if there is no mass loss. 
Present observations have observed clumps more massive than 
this criterion in clumpy galaxies. Therefore the existence of such 
massive clumps is not unrealistic. However, current high red- 
shift observations are biased to relatively luminous galaxies, and 
such brighter galaxies are naively expe cted to host more mas- 
sive c lumps than fainter galaxies (e.g.. lElmegreen et alJ 12007. 
2009a,3). Hence, although the massive clumps can exist, their 
typical mass may be lower than the critical mass, and the clump 
masses could depend on the mass of the host galaxy. 

In 33.41 I showed that the life time of clumps is also a 
key factor. Although I arbitrarily set the mass loss rate of the 
clumps (99 % mass los s in 0.3 5 Gyr) in my calculation, sim- 
ulations of iGenel et all (120121) have suggested much shorter 
life times of clumps, ~ 50 Myr at R ~ 3 kpc. Therefore, I 
infer that the kinematic imprint would be even weaker than 
my calculation or could not be seen if such a high mass loss 
rate is real. However, it is still being actively debated whether 
the clumps are short or long lived, as I noted in 93.41 (see 
Genzel et alj|201 U iNewman et alj|2012l: iForster Schreibe r et alJ 
201 HlGuo et al.H2012l) . 

Other dynamical processes may be able to provide a halo 
system with angu l ar m omentum. Numerical simulations by 
iBekki et all d2002l 120051) have demonstrated that metal-poor 
GCs can obtain a net rotation by galactic mergers. However, 
their results showed that the GC system rotates significantly in 
an outer halo region rather than in an inner region. Moreover, 
the direction of the rotation would not n ecessarily coincide with 
the disk rotation in the merger scenario. iBois et all d201 lh con- 
ducted binary merger simulations with large samples, which also 
indicate that merger remnants generally have larger spin param- 
eters in their outer regions t han in their inner regions. In addi- 
tion, iDi Matteo et alJ (12009) demonstrated that a merger rem- 
nant of non-rotating spherical galaxies can be a rotating ellip- 
tical. However, also in their simulation, an outer region of the 
remnant rather than an inner region, indicates significant rota- 
tion. In Figure [9] I show ratios of LoS velocities to LoS veloc- 
ity dispersions. The ratios decline sharply outside the disk radii, 
and the outer regions are not rotating. Thus the rotation signa- 
tures by clumpy disk formation are different fro m those induced 
by gal actic mergers (compare to Figure 6 of IDi Matteo et al] 
(2009)). Also in cosmological A^-body simulations, DM halos 
generally have averaged rotations approximately proportional to 
j oc r u (e.g., Bu llock et alj|200lb . This indicates nearly con- 
stant or slightly increasing rotation velocities toward outer halos, 
where j is specific angular momentum, although a stellar halo 
may not necessarily share the same kinematic state with a DM 
halo. In FigureQX)] radial profiles of angular momenta in my cal- 
culations are shown. Obviously, the angular momentum profiles 
are not in monotonic fashion and show little rotation outside the 
disk radii. I therefore suggest that the disky halo rotation caused 
by clump motions could be unique and distinguishable from the 
cosmological merger-origin rotation. 

Having formed in a galactic disk, old stars may also have 
properties similar to the HSs influenced by clumps, even if the 
disk has not experienced a clumpy phase. These old stars can 
be heated up in the vertical direction by perturbations such as 
passages of spiral arms and giant molecular clouds. Eventually 
the old disk stars can have a thickened density distribution. 
However, the disky halo rotation by clumps is much slower than 
disk rotation. As a result, the rotating HSs could be distinguished 
kinematically (and also chemically) from the old disk stars. 

Meanwhile, GCs associated with a galactic disk (disk 
GCs), which have been observed in some galaxies, also have 
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Fig. 9. Profiles of mean LoS velocities over dispersions of the 
halo systems along the disk planes. The LoS dispersion is com- 
puted in each bin, not in the central region. The red solid line 
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10. Spherically averaged specific angular momentum pro- 
of the halo systems. The line types are the same as in Figure 



a disky kinematic state similar to my result (jZinnl f l985: 
I Chandar et al.l l2002t iMorrison et all l2004t iMackev & Gilmord 
l2004t iMackev & van den Berghll2005l) . However, the disk GCs 
are metal-rich and thought to have formed in a relatively late 
stage of galaxy evolution, although there seems to be a hand- 
ful of old and metal-poor disk G Cs (e.g.. iDinescu et al 1 I1999L 
120031: IMackev & Gilmord l200H) . Furthermore, their orbits al- 
most coincide with disk rotations, indicating nearly circular ro- 
tation velocities and low orbital eccentricities, which are incon- 
sistent with my results (Figure |4] slow rotation and mild orbital 
circularization). Thus, I expect that the disk GCs are not indica- 
tive of the clump interaction and that clumpy disk formation does 
not explain the origin of the disk GCs. 

However, origins of old GCs and HSs have not been 
elucidated yet and their initial st ates of distribution an d 
kinematics are thus unclear (see Brodie & Strader 2006). 
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Some studies have suggested that a part of GCs and HSs 
seem s to have external origins from accreting dwarf galaxie s 
(e.g.. ICote etaLlll998h lAbadi et all 12005 iMackev et al1l2010b . 
Recently, iTrujillo & Bakosl (1201 2h observed stellar halos of two 
Milky Way-like galaxies at a redshift of z ~ 1, which seem 
to have structural properties similar to those of counterparts in 
the local universe. They suggest, however, that some accretion 
and/or star formation is still lasting after z ~ 1 . If it is the case, 
acquisition by a galaxy of GCs and HSs may postdate the clumpy 
disk formation and the halo objects originating from such late 
accretion can avoid the influence in the clumpy disk formation 
stage. 

What fraction of disk galaxies can have the imprints of past 
clumpy phases? If it is here assumed that a fraction, A, of current 
disk galaxies have experienced clumpy disk formation and that 
a fraction, B, of the clumpy galaxies have hosted giant clumps 
of M c i > 5 x 10 8 M , it is simply expected that approximately 
A x B of the current disk galaxies would show the imprinted 
signatures, even though the clumps have already disappeared to 
the current galaxies. However, it is still challenging to estimate 
the fractions of A and B in present observations that would be 
biased toward luminous galaxies. In this regard, cosmological 
simulations may be powerful tools to explore this issue, although 
curre nt simulations are s till lacking samples with high resolution 
(e.g.. lMartigeTaT1l2012l) . 

Are there remnants of clumpy galaxies nearby? Some local 
disk galaxies have indeed been observed to have a net rotation in 
their old (m etal-poor) GC system , such as th e Large Magellani c 
Cloud (e. g iFreeman et al.1 1983b . M31 (e. gJPerrett et ajj2002|) . 
NGC253 dOlsen et all 12004. N GC524 dBeaslev et al.1 120041) 
NGC 2683 dProctor et al.1 120081) . NGC3115 dKuntschner eTall 
2002)Q However, because of their faintness, the number of 
GCs observed spectroscopically to measure their kinematics is 
still limited in extra-galaxies. Furthermore, it would be obser- 
vationally difficult to pick out a sufficiently large number of 
GCs from their disk regions. In the Milky Way, there seems 
to be little net rotation i n me tal-poor GCs in an y radi a l rang es 
(e.g., iBrodie & Straderl d2006lFl although see izinnl dl985lf ). 
Meanwhile, rotation of HSs in the Milky Way is still contro- 
versial. Some observations have suggeste d a slow prograde ro- 
tation, < 50 kmsT 1 , of the inner halo (rChiba & Beerslfe OOO ; 
iKeolev et al.ll2007tlCarollo et al.ll2 007: Sna gna et al.ll2010H 5 L but 
others have favored a non-rotating inner halo ( Val lenari et al.l 
2006: ISmith et al.ll2009t [Bond et alJl2010l) . Thus, present obser- 
vations do not seem to enable us to identify galaxies that were 
once clumpy disks. A near-future astrometry satellite mission, 
Gaia, is promising in this regard. 

It should be noted that my calculations were performed with 
the simplistic toy models. In a real galaxy, disk formation and 
growth of DM halo follow the clumpy phase. A massive disk 
can deform a galactic potential into an axisymmetric (disky) one. 
Then, a distribution of halo objects would become aspherical, 
like Figure [2] by the disk potential even if there were no clumps 
in a disk-formation s tage. M oreover, using a numerical simula- 
tion llnoue & Saitoh! (1201 ll) have demonstrated that a DM halo 
initially assumed to be an NFW profile is turned into a cored 



3 In M33, GCs older than 1 Gyr have been observed to have a sig- 
nificant rotation, but the nu mber of observed sam ples older than several 
gigayears is only a handful (Chandar et al. 2002). 

4 However, Brodie & Stradei] d2006h have mentioned that a strong 
prograde rotation is seen in the most metal-poor GCs. 

5 ICarollo et al. ( 2007) did not reject the possibility of a non-rotating 
inner halo. 



profile in a clumpy galaxy, though I assumed the rigid NFW halo 
potential. More sophisticated and realistic simulations should be 
performed in future works. 

5. Summary 

If halo objects such as GCs and HSs have formed by internal 
processes of galaxy evolution, clumpy disk formation can be ex- 
pected to have taken place in a pre-existing halo consisting of 
the old objects. In such a case, the clump orbital motions may 
be able to exercise influence on the kinematic state of the halo 
objects. Accordingly, I performed restricted A^-body calculations 
and studied the kinematic imprint of the clumpy phase on a halo. 
I found that the halo objects can be expected to indicate slow ro- 
tation and orbital circularization around the galactic disk after 
the clumpy phase. I also found that postulated clump masses 
must be more massive than 5 x 10 8 M and five clumps are 
enough to leave the imprint. 
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